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Motivated by the recent discovery of high temperature antiferromagnet SrRu206 [1, 2] and its 
potential to be the parent of a new superconductor upon doping, we construct a minimal t 2 g- 
orbital model on a honeycomb lattice to simulate its low energy band structure. Local Coulomb 
interaction is taken into account through both random phase approximation and mean field theory. 
Experimentally observed antiferromagnetic order is obtained in both approximations. In addition, 
our theory predicts that the magnetic moments on three t 29 -orbitals are non-collinear as a result of 
the strong spin-orbit coupling of Ru atoms. 


I. INTRODUCTION 

Magnetism and superconductivity are closely related 
to each other, as a common thread in several families 
of unconventional superconductors. [3] Singlet Cooper 
pairing is expected to be mediated by antiferromagnetic 
(AF) fluctuations near the AF phase boundary like in 
most cuprates [4], iron-based [5], and heavy fermion [6] 
superconductors. While triplet Cooper pairing is widely 
believed to be triggered by ferromagnetic fluctuations 
as in Sr 2 Ru 04 [7-9]. Therefore, looking for unconven¬ 
tional superconductivity in materials with strong mag¬ 
netic fluctuations is a guiding principle in the commu¬ 
nity. SrRu 206 is synthesized recently [1] and reported as 
an antiferromagnet with a Neel temperature as high as 
565K [2]. Similar to other ruthenate Sri+„Ru„Oi+ 3 „ [10], 
SrRu 2 06 is also a layered material. Due to the layered 
property of the sample, the AF order may be easily de¬ 
stroyed by introducing quantum fluctuations via doping 
or high pressure. As a result, SrRu206 may be a good 
parent compound to realize high temperature supercon¬ 
ductivity. 

Different from ruthenate Sri+„Ru„Oi-|- 3 „ in which the 
RuOg octahedra are point-sharing and form a square 
lattice [11], in SrRu 206 the RuOe octahedra are edge¬ 
sharing and the Ru atoms are arranged on a honeycomb 
lattice. According to first-principle calculations, [2, 12] 
the t 2 g"Orbitals of Ru are found to dominate the low en¬ 
ergy states. This suggests that a t 2 g-orbital model on a 
honeycomb lattice is a relevant minimal model for fur¬ 
ther studies on the correlation effect. On the other hand, 
even though the spin and charge degrees of freedom have 
been broadly studied in the honeycomb lattice (such as 
graphene) with only 7r-electrons [13], the orbital degrees 
of freedom would bring us new features. For instance, 
the studies of (pa;,Py)-orbital models on the honeycomb 
lattice revealed Wigner crystallization [14] and anoma¬ 
lous quantum Hall effect [15, 16]. In parallel, SrRu206 
provides us a natural realization of the Og d-orbital on 
the honeycomb lattice. 

In this paper, we first derive an effective t 2 g-orbital 
tight-binding Hamiltonian as a minimal model for 
SrRu 2 06. We then consider the correlation effect 
through both random phase approximation (RPA) and 



FIG. 1. Schematic discription of the lattice structure and all 
hopping elements up to the 3rd nearest neighbour. Three t 2 g- 
orbitals of Ru are placed on the honeycomb lattice. O atoms 
(small circles) are distributed above (yellow) or below (red) 
the Ru-plane. Each RuOe octahedron is slightly distorted. 

mean field theory. We obtain the experimentally ob¬ 
served AF order and estimate the Neel moment and tran¬ 
sition temperature within the mean field theory. Further¬ 
more, we find the orbital-resolved AF moments on three 
t 2 g-orbitals are non-collinear as a result of the strong 
spin-orbit coupling (SOC) on Ru atoms. Our minimal 
model provides the basis for further theoretical studies, 
and the magnetic structure we uncovered would trigger 
further experimental interests in this new member of the 
ruthenate family. 


II. t 29 -ORBITAL TIGHT-BINDING 
HAMILTONIAN 

In SrRu206, the low energy bands mainly come from 
the t 2 g-orbitals of Ru atoms as reported by first-principle 
calculations. [2, 12] The Ru atoms form a honeycomb 
lattice in each RuOe layer. Therefore, to mimic the sys¬ 
tem by a minimal model, we only consider the t 2 g-orbital 
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electrons on a honeycomb lattice as shown in Fig. 1. The 
coordinate is set up with the origin at Ru-site and three 
axes point to three oxygen atoms above the Ru-plane in 
the undistorted RuOe octahedron. In this coordinate, 
the c-axis perpendicular to Ru-plane is along (1,1,1)- 
direction, and the dxy, dxz and dyz orbitals forms an ex¬ 
act t^g multiplet, as schematically represented in Fig. 1. 
The blue/yellow/red lobes lie within/above/below the 
Ru-plane. Although in SrRu 206 , the RuOe octahedron 
is slightly twisted around and stretched along the c-axis, 
the t 2 g-orbital degeneracy is protected by the sym¬ 
metry. As a result, the coordinate established in the 
undistorted octahedron will be used in this work. 

To construct a tight-binding Hamiltonian, we keep 
nearest neighbor hopping next nearest neighbor 

hopping t 5 ,-., 8 , 5 ', 8 'i and the third-neighbor hopping tyr^ 4 r. 
All these hopping elements are schematically shown in 
Fig. 1. We notice that 7 ^ tg and ts 7 ^ *iae 
to the distortion of the RuOe octahedra. Further¬ 
more, we add intra- and inter-orbital on-site energies 
Vo and bg', respectively. In a pure two dimensional 


model, Vq should be exactly zero due to the orthog¬ 
onality of different orbitals. But here we are looking 
for an effective two dimensional model in the kc = 0 
plane and thus the inter-layer hopping could lead to an 
effective on-site inter-orbital mixing Vq. In addition, 
we add an SOC term Hsoc = ® (Jy.i’i, 

where ' 0 ^ — [di^xy,^^di^xz,^jdi^yz^-^^di^xy^^Tdi^xz,].^di^yz^^\^ 

o'fi = [CTa;,cTj,,cr^] are three Pauli’s matrices, and = 
[Lxi Ly, Lz] are rank-3 angular momentum matrices act¬ 
ing on orbital space, with all nonzero elements given by 

Lx,12 = —Lx,21 = —*) Ly^io = —Ly,31 = —i, 

Lz,23 = —Lz,32 = L (1) 

With the above ingredients and for a given inter¬ 
layer momentum kc = 0 , we arrive at an effec¬ 

tive two-dimensional tight-binding Hamiltonian H = 
Lk ^k^k^'k in the basis = [V'Akt, V’skt, ' 0 Ak;, V’sk;] 
where V’*,ka = [4,xy,ko-, ds,xz,i^a, ds,yz,\^a] for s = A, B on 
the two sublattices. The matrix is explicitly written 
as 
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where (ai, a 2 , ag) are three displacements from an A-site 
to its nearest neighbour B-sites. All the model param¬ 
eters are determined by fitting the first-principle band 
structure [12]: (in unit of eV) ti = 0.16, ^2 = —0.01, 


I 

ts = 0.30, t4 = -0.02, to = -0.10, t'o = 0, to = -0.01, 
tr = 0.04, ts = 0.11, t's = 0.02, t '4 = -0.04, = -0.01, 

-0.01, = -0.01, Vo = -0.09, Vg' = -0.07, 

0.16. In particular, we obtain the SOC strength 
0.16eV, which is in agreement with the literature 
[17]. Based on these parameters, the band structure of 
our minimal model is plotted in Fig. 2. 

In the paramagnetic state SrRu 2 06 is a band insula¬ 
tor, as seen in our band structure and the first-principle 
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FIG. 2. Band dispersion of our f 29 -orbital tight-binding model 
along high symmetry lines. The local-density-approximation 
result (from Ref. 12) is also shown (red symbols) for compar¬ 
ison. 



FIG. 3. (a)Vertex function T(afj,n;biym) defined in Eq. 6 in 
the magnetic channel. (b)The Feynman diagrammatic repre¬ 
sentation of the magnetic susceptibility matrix (Eq. 7) within 
the RPA. 


result [2, 12]. However, a strong Hund’s coupling would 
bind up the electrons to form a spin-3/2 state. Due to 
the bipartite lattice an AF order with moment d/rs/Ru 
is expected. Such an AF order has already been observed 
by Neutron scattering experiment in SrRu206,[2] but the 
observed moment is only l.S/is/Ru, much smaller than 
SfiB- This indicates the inadequacy of a naive local mo¬ 
ment picture. Instead, the itinerant property of electrons 
and SOC may play important roles. In the following, we 
will investigate the effect of correlation and SOC on the 
AF order in the itinerant picture of the t 2 £("Orbital model. 


III. INTERACTION AND 
ANTIFERROMAGNETISM 

We adopt general multi-orbital local Coulomb interac¬ 
tions as follows, 

i a a>b 

+ J (^lu^io-blg./CLia' + J ] ) (5) 

a>b,(T(T' a^b 

where rua = u is the Hubbard 

interaction, V is the inter-orbital charge interaction, J 
is the Hund’s coupling and J' is the pair hopping term. 
These four interactions satisfies the relation J' = J and 
U = V + 2 J. [18] Among these four terms, only U and 
J are responsible for magnetic channel instabilities. [19] 
So in the following discussions, we will only retain the U 
and J terms. 

Since the non-interacting model is a band insulator, 
the bare susceptibility only depends on momentum very 
weakly. So we perform an RPA level calculation instead, 
since RPA will pick out relevant channels and strongly 
enhance their susceptibilities. 


The interactions {U and J terms) are first written in 
magnetic channels: —SiT(I;J)Sj, where I = (afin) and 
J = (bvm) with the orbital index {a/b), spin-direction 
index {fJ.ji'), and sublattice index (n/m). The spin op¬ 
erator is = 5 Sa /3 «laCra/3«n/3) ^nd the vertex 

function T[a^n]bvm) represented by Fig. 3(a) is diago¬ 
nal in both spin-direction and sublattice subspaces and 
is given by 


f (a/in; bvin) 


nm 1 ^ ^ 

d ^ b 


( 6 ) 


The RPA is a bubble summation as represented in 
Fig. 3(b). After solving the iterate equation we obtain 
the magnetic susceptibility matrix as 


X(q,w„)= i-fx(°)(q,in„) 


x(°^(q,zn„). (7) 


,ivn) is the bare susceptibility whose matrix ele¬ 
ment is defined as 

j-l/T 

Xui^,i^n)= / (S'/(-q,r)S'j(q,0))e*'"’'^dr. 

( 8 ) 

Here, we have used Matsubara frequency with T the tem¬ 
perature. 

We use the interaction parameters U = 1.35eV and 
J = 0.14eV, being half of the values obtained by first 
principle calculation [2], since RPA is known to over¬ 
estimate the magnetic instability. We plot the leading 
(largest) eigenvalue A of the hermitian susceptibility ma¬ 
trix x(q, 0) as a function of q in Fig. 4. The peak at 
q = 0 implies a magnetic instability that is periodic 
across the unit cell. The corresponding leading eigenvec¬ 
tor decides the form factor of the magnetic order, namely 
the magnetic structure within a unit cell. This is shown 
in Fig. 5(a). We find the total moment (of the three or¬ 
bitals) is along the c-axis and changes sign from one to 
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FIG. 4. The largest eigenvalue A of = 0) along high 

symmetry cnts. Here U — 1.35eV, J = 0.14eV, ksT = 0.3eV. 


the other sublattice within the unit cell. This is exactly 
the AF order seen in the neutron scattering experiment. 
More interestingly, the moments on the three orbitals are 
non-colinear about the c-axis. This is a prediction of the 
present work. 

Next, we employ mean field theory to quantitatively 
investigate the AF moment and the transition tempera¬ 
ture. The interactions are decoupled in magnetic chan¬ 
nels: 

41/ 2/7 

^-^(S.a)-S,, + —(S,,).(S,„), (9) 

and 

J ^ ^ ' ^ib 

gg' 

^ -2J(S,,) • - 2JS,, • (S,b) + 2J(S,,) • (S,t,). 

( 10 ) 


Then we obtain the mean field Hamiltonian 




-2 J ^ (S,,) • S, 

i,a^b 


( 11 ) 


where Hq is the non-interacting part, which is the same 
as Eq. 2 but written in real space. The order parameters 
(Sia) = I Z)a/3(®L'^a/3«»/3) ^re then numerically solved 
iteratively until convergence is achieved. 

Our mean field result confirms the non-collinear AF 
configuration [Fig. 5(a)] revealed by RPA in the normal 
state. We have performed unrestricted mean field calcu¬ 
lations starting from random initial spin configurations. 
No translation symmetry is assumed in advance. How¬ 
ever, the results all converge to the same non-collinear 
AF configuration up to a shift of the sublattice. The 


(a) total moment 




FIG. 5. (a)The AF configuration is schematically shown 
within a unit cell, with three orbital resolved moment rep¬ 
resented by three different colors. The total moment is along 
c-axis [(1,1, l)-direction]. (b)The total AF moment M and 
the deviation angle 9 vs temperature within the mean field 
approximation. 


magnetic unit cell is always found to be equal to the 
lattice unit cell, which is in agreement with the unique 
peak q = 0 in the leading eigenvalue of the RPA mag¬ 
netic susceptibility in the momentum space as shown in 
Fig. 4. The non-collinear AF configuration is the result 
of the strong SOC on Ru atoms. The values of the to¬ 
tal moment M and the deviation angle 9 of each orbital 
moment relative to c-axis [shown in Fig. 5(a)] vs temper¬ 
ature are plotted in Fig. 5(b). Two sets of parameters 
are used. One is (C/, J) = (2.7, 0.28)eV from first princi¬ 
ple calculation [2] , and the other is weakened by a factor 
of two, (C/, J) = (1.35,0.14)eV. The latter leads to re¬ 
duced M and transition temperature. Since mean field 
theory overestimates the ordering, further considerations 
of quantum/thermal fluctuations beyond the mean field 
theory are necessary for quantitative comparison to ex¬ 
periments. Interestingly, we find 9 becomes smaller with 
increasing interaction strength. This is because a larger 
J tends to align the spins on different orbitals, while 
a larger SOC breaks the Hund’s rule more significantly. 
Since neutron scattering only ’sees’ the total moment, we 
expect this particular kind of orbital-resolved AF order 
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can be observed in more delicate experiments like orbital- 
selective nuclear magnetic resonance[20, 21] through the 
anisotropic hyperfine interaction [22] or orbital-resolved 
angle-resolved photoemission spectroscopy[23] through 
the photon polarization selection rule [24]. 

IV. SUMMARY AND FUTURE WORKS 

In summary, we have constructed a t 2 g-orbital model 
on a honeycomb lattice. Local Coulomb interaction was 
investigated in both RPA and mean field theory. Ex¬ 
perimentally observed Neel order is obtained. Further¬ 
more, our theory predicts that the magnetic moments on 
three orbitals are non-collinear as a result of the strong 
spin-orbit coupling of Ru atoms. This particular kind of 
orbital-resolved AF order is expected to be observed in 
future experiments. 

For future works, possible superconductivity in this 


compound after doping or under pressure is an interest¬ 
ing direction. Our Og-orbital model can be used as a 
minimal model to study possible superconductivity upon 
doping. The AF fluctuation may induce singlet Cooper 
pairing between the nearest neighbours. However, due 
to the strong SOC, triplet pairing may coexist with the 
singlet pairing. This part of work is being in progress. 
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